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Abstract 

To elucidate the role played by the transition metal ion in the pnictide materials, we compare 
the electronic and magnetic properties of BaFe2As2 with BaMn2As2. To this end we employ the 
LDA+Gutzwillcr method to analyze the mass rcnormalizations and the size of the ordered magnetic 
moment of the two systems. We study a model that contains all five transition metal 3d orbitals 
together with the Ba-5d and As-4p states (ddp-model) and compare these results with a downfolded 
model that consists of Fe/Mn d-states only (d- model). Electronic correlations are treated using the 
multiband Gutzwiller approximation. The paramagnetic phase has also been investigated using 
LDA+Gutzwillcr method with electron density self-consistency. The renormalization factors for 
the correlated Mn 3d orbitals in the paramagnetic phase of BaMn2As2 are shown to be generally 
smaller than those of BaFe2As2, which indicates that BaMn2As2 has stronger electron correlation 
effect than BaFe2As2. The screening effect of the main As 4p electrons to the correlated Fe/Mn 3d 
electrons is evident by the systematic shift of the results to larger Hund's rule coupling J side from 
the ddp-model compared with those from the d-model. A gradual transition from paramagnetic 
state to the antiferromagnetic ground state with increasing J is obtained for the models of BaFe2As2 
which has a small experimental magnetic moment; while a rather sharp jump occurs for the models 
of BaMn2As2, which has a large experimental magnetic moment. The key difference between the 
two systems is shown to be the d- level occupation. BaMn2As2, with approximately five d-electrons 
per Mn atom, is for same values of the electron correlations closer to the transition to a Mott 
insulating state than BaFe2As2. Here an orbitally selective Mott transition, required for a system 
with close to six electrons only occurs at significantly larger values for the Coulomb interactions. 
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I. INTRODUCTION 



The discovery of high Tc superconductivity in LaOi-j^F^^FeAs has initiated a detailed 
investigation of these and related iron based compounds with the ultimate aim to identify 
new high temperature superconductors Many binary, ternary and quaternary compounds 
have been synthesized and investigated experimentally, some of which were quickly found 
to exhibit high T^ superconductivity with electron or hole doping or under pressure or 
intrinsically |2|-l7| • Meanwhile, several closely related and physically interesting systems have 
also been studied. An interesting example is the BaMn2As2 compound, which has the same 



layered tetragonal ThCr2Si2-type crystal structure as the BaFe2As2 at room temperature 
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BaMn2As2 is an antiferromagnetic (AFM) insulator with a small gap; while BaFe2As2, 
a parent compound of the Fe-based high Tc superconductors, is an AFM metal at low 
temperatures. The more localized behavior of the ordered magnetic state in the Mn-based 
system suggests that electronic correlations in this system are stronger or more efficient in 
changing the electronic properties of the material. A closer, comparative investigation of 
distinct transition metal ions therefore promises to reveal further whether iron is indeed 
special among its neighbors in the periodic table. This is also expected on quahtative 
grounds on the basis of the Hund's rule picture [l3|. 

Theoretically, the description of the normal states for iron pnictide superconductors is 
complicated by the multi-orbital nature and the presence of electron correlation effects. 
Interesting results can be obtained from low energy theories that are primarily concerned 
with states in the immediate vicinity of the Fermi surface 14-19|. However, material specific 



insight often requires a careful analysis of states over a larger energy window. In the iron 
pinctide and chalcogenides the strength of the correlations is controlled by the Hund's rule 
coupling J, rather than by the Hubbard U or the p-d charge transfer energy 13j. This 
makes these materials a new class of correlated metals, different from doped Mott or charge 
transfer insulators. Calculations based on density functional theory (DFT) usually yield too 
wide Fe 3d bands and too large magnetic moment in the AFM ground state 12||. Recently 
the AFM phase of the iron pnictide parent compounds has been successfully described by 
using a combination of DFT and dynamical mean field theory (DMFT) with full charge 



self-consistency 



201 ]. These investigations clearly demonstrated the important role played by 



the Hund's rule coupling J for the magnetic and electronic properties of the iron pnictides. 
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A systematic analysis of a large class of Fe-based materials further demonstrated that the 
nature and strength of local correlations are rather universal, while material specifics result 
from changes in the kinetic energy, amplified by^ the frustrated nature of electron hopping 
path due to Fe-As-Fe and direct Fe-Fe overlaps[21]. This universality in the value of the local 
Coulomb interactions U, U', J only exists if one includes in addition to the Fe-3d states at 
least the As 4p orbitals. Downfolding to models with iron states only is sufficient to lead to 
material specific variations in the local correlations. In view of these insights is it important 
to compare the behavior of Fe-based systems with systems based on other transition metals, 
also in view of the ongoing debate whether these systems should be considered weakly or 
strongly correlated. 

On a technical level, quite a few efforts have been devoted to study the electronic and 
magnetic properties of the systems using the Gutzwiller variational method {22!] and simpli- 



fied tight-binding models 23|, |2J| with the particular appeal that the Gutzwiller approach 



is computationally cheaper if compared to more sophisticated DMFT approaches such as 
the usage of Quantum Monte Carlo algorithms. For example, the experimental magnetic 
moment was obtained in a three band Hubbard model of LaOFeAs which contains only Fe 



3d t2g orbitals 23|]. Calculations based on the five-band model show that the Cg orbitals are 
also important and the system exhibits a sharp transition from paramagnetic (PM) state to 
AFM state with large magnetic moment when J increases, which is different from the results 
of three-band model [2^. Another important question is whether, in addition to the Fe 3d 
Cg orbitals, the inclusion of As 4p elections is required for a more quantitative description 
of the AFM and paramagnetic phases. Finally, recent calculations for LaOFeAs based on 
the LDA+Gutzwiller approach that includes a Fe3d-As4p Wannier-orbital basis does indeed 
reproduce the experiment ally observed small ordered magnetic moment over a large region 
of (f/, J) parameter space |25[|. 

In this paper, we report a comparative study of the BaFe2As2 and BaMn2As2 on the 
electronic and magnetic properties by investigating the Fe/Mn d-band model and the ddp- 
model that includes Fe/Mn 3d, Ba 5d, and As 4p orbitals with Gutzwiller approximation 26- 



28| which has been shown to be equivalent to a slave-boson mean field theory 29l|. The 



paramagnetic phase for the two systems has a. 
method with electron density self-consistency 



so been investigated using LDA+Gutzwiller 
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3l|. BaMn2As2 is shown to be a more 



correlated material than BaFe2As2. The screening effect of the As 4p electrons is evident 
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by the comparison between the results of the Fe/Mn d-model and those of the ddp-model 
or LDA+Gutzwiller method. The key difference in the behavior of these systems is shown 
to be the distinct electron count of the 3d states, a result that is robust against changes 
in the detailed description of these systems. While the even number of electrons for an 
Fe site requires that the system undergoes an orbitally selective Mott transition, Mn with 
five electrons can more easily localize and therefore has a much stronger tendency towards 
large moment magnetism and Mott localization, even for the nominally same values of the 
Coulomb interactions. The prospect to tune the tendency towards localization by varying 
the transition metal ions may therefore play an important role in optimizing the strength of 
electronic correlations for high temperature superconductivity. 



II. FORMALISM 



For the sake of being self-contained, we outline the main formalisms for the multiband 
model with the Gutzwiller approximation 26N28| and the LDA+Gutzwiller approach with 



electron density self-consistency |30l. |31|. 



A. Multiband model with the Gutzwiller approximation 

The multi-band Hamiltonian with local many-body interactions is written as 

•H = -Ho + 'Hi (1) 
where the bare paramagnetic band Hamiltonian (7/o) is 

T^O = ^ ] ^iaj pC^^^Cj j3a + ^ ^ ^iaCiaa^iaa (2) 

tiaj/3 is the electron hopping element between orbital at site i and orbital ipp at site j. Sia 
is the orbital level. For the Fe/Mn d-model, a and /3 run over the Fe/Mn 3d orbitals. For 
the ddp-model, a and /3 run over the Fe/Mn 3d orbitals, Ba 5d orbitals, and As 4p orbitals. 
c\c) is the electron creation (annihilation) operator, a is the spin index. 
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The typical Hubbard density- density type interaction term (Tii) for Fe 3d electrons is 



'Hi = f/ ^ Ui^^rii^l + U' ^ 

-J E 

2,7<7',(T 



Tli-yQ-Tli'yf (jf 



(3) 



which is usually sufficient for describing collinear magnetic order 
Fe/Mn 3d orbital index, U' = U — 2J, and riiya = c^i^a'^i-fa- 
We take a variational wave function of the Gutzwiller form, 

Gl^o) 



23| . Here 7 is the correlated 



^2 



(4) 



with the Gutzwiller approximation to calculate the expectation values of Eq. [T] 26|-|28| . ^'o 
is the uncorrelated wave function. G is the Gutzwiller projection operator, which adopts 
the following form 

G = e-^'^^»^l-^'><^»l (5) 
where is the Fock state generated by a set of |c|^o-| 



(6) 



with n: 



T 



or 1, which identifies whether there is an electron with spin a occupied in 
orbital 7 for Fock state T at the z*'* site, gu is a variational parameter which controls the 
weight of the local Fock state | J-^) in the Gutzwiller wave function, = 1 for the empty and 
singly occupied configurations because in these cases there are no electron-electron repulsion 
involved. According to Ref. 28|, the expectation value of the electron Hamiltonian "H can 
be expressed as 



ia,jl3,cr 



where lc\cjR^ 



^0 



orbital renormalization factor Zi, 



(7) 



\l/oy and Uijr = {J^i\l-Li\ J^i) . We define the Gutzwiller 
E 1 for Ba 5d and As 4p orbitals in the ddp-model since 
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they are not considered as strongly correlated orbitals. Following Ref. 28|, for the Fe/Mn 
3d orbitals 



1 



770" 



(8) 



with n\^^ 



{ni-ya)^- ViT is the occupation probability of configuration which absorbs 
Qij: and serves as the additional variational parameter due to the Gutzwiller wave function. 

Since in the Hamiltonian of Eq. [T]is obtained by downfolding the DFT band structure, 
some contributions from "Hi have already been taken into account by "Hq in a mean-field way. 
Such contributions are commonly referred to as the double counting term. Consequently, 
the total energy of the system is given by 



with 



Etot = {^)g + 



Ed.c. = -Yl {UN^d {Nid -1)12- JN^d {N,d/2 - 1) /2) 



(9) 



(10) 
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following the treatment in DFT+U method 
(Hund's rule coupling J) parameter, which for d orbitals is given as 



33|. Here U (J) is the averaged Coulomb U 
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34| 



(11) 
(12) 



j = U -U' + J 

^id = Yli'ya''^%c7 expectation value of the total Fe/Mn d-orbital occupancy at site i. 

The double counting term E^.c. introduces some spin-independent chemical potential shift 
for the Fe/Mn 3d orbitals. 
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-U ( iV,, - 1") + J (n,J2 - i 



(13) 



In the Fe/Mn d-model that contains only 3d states, Si^ has no effect since it merely shifts 
all the Fe/Mn 3d levels by a constant. However, in the ddp-models e-i^ becomes important 
because it changes the position of the Fe/Mn 3d levels with respect to Ba 4d and As 4p 
orbitals and hence affects the electron flow between these orbitals. 



The minimization of Etot with the constraints of YlirPi^ ~ 1 ^"^^ X]j-P«^'^i7cr 
yields an effective single-particle Hamiltonian 



"'t-ycr 



(14) 



where fii 

and /ijacr 



dzj 



+ rjiaa if a belongs to the correlated Fe/Mn 3d orbitals 



= Eia otherwise. It has been shown that Hes describes the Landau-Gutzwiller 
quasiparticle bands 37[ , and the square of the orbital renormalization factor z"^ corresponds 
to the quasiparticle weight [31j. {zi^„} and {r^j^o-} are obtained by solving the following 



equations, 
where 



70" 















n) 





(15) 



(16) 



with = Y^jp yzjiB^Uajfs {cl^^Cjp^j ^ + c.c.j. \l'o froui Efl,^ [Hand {pir} from Eq. [I5]need 
to be solved self-consistently to obtain the final solution 35 1. 



B. LDA+Gutzwiller 



The above formalism based on the Gutzwiller approximation on the multiband model 



can be naturall y c ombined with DFT, in which full electron density convergence can be 



achieved 



3q|,|3i[ 



The total energy density functional can be written as 



E[p] = {^Gf^G) + EH [p] + ^.c + / P (r) v,,n (r) + E, 



-'ion— ion 



(17) 



where T is the kinetic energy operator and Eh is the Hartree potential energy. The electron- 
ion and ion- ion potential energies are given by / p (r) Vion (r) d^r and Eion-ion, respectively. 



In the LDA+Gutzwiller method 
takes the following form 
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3l|, the approximate exchange-correlation energy E^ 



E,, = E^^""^ [p] + Y U^TV^T + ^d.c. 



(18) 



The electron density p under the Gutzwiller approximation can be written as 



P (r) = Y ^^c.a^^j/^aV'L (r) Vjp (r) (c^Cj/j^)^ + Y l'^^" {niaa)o (19) 

ia,jl3,a 



where {ipia} is a general complete basis set including local correlated orbitals (e.g., Fe/Mn 
3d orbitals) as well as nonlocal orbitals. indicates that the summation excludes the term 
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z X r z 



FIG. 1: (Color online) Band structure of BaFe2As2 in PM phase based on LDA (solid line), Fe 
d- model (circle) and ddp- model calculations (cross). Fermi level is shifted to zero. 

with ia = jp. The kinetic energy is also renormalized 



G 



T 







The total energy E [p] of the closed form can be similarly minimized as in the above multi- 
band model with respect to \&o and {pij-}. The resultant effective single-particle Hamiltonian 
has the same form as Eq. [T3] with the hopping element 

Ua,p = (vrc \f\ ¥^,73) + j (r) + v,on (r)) </^L (r) (r) d\ (21) 

The LDA screening potential v^^^ (r) is the sum of the Hartree potential vh (r) = ^^f^ 

and the exchange-correlation potential f^c^^ (r) = ^^^^f^-^- The renormalized orbital level 

A^iacr = ^ia + ^ia + Ctacr + ^iau if oc bclougs to the Correlated Fe/Mn 3d orbitals and 

( / t \ \ 

^^^acT = ^ia OtherwisC. Here £ia = tiaia and Cj^o- = \ y^i-^icT^jP'^l ^ ~^ C.C.j . {Ziya-} 

and {rjiya} are obtained by solving the same equation as Eq{T5l One can see that the major 
difference is that the hopping element tiajp and the orbital level will be updated self- 
consistently by the LDA screening potential v^^"^ (r) in LDA+Gutzwiller method, while 
they are fixed in the multiband model. 

III. EFFECTS CAUSED BY DOWNFOLDING 



The Fe/Mn d-model and ddp-model for BaFe2As2 and BaMn2As2 have been obtained by 
downfolding the first-principles LDA electron bands of the experimentally determined crystal 



r z 



FIG. 2: (Color online) Band structure of BaMn2As2 with settings same as Fig. [TJ 
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FIG. 3: (Color online) The atomic configuration of the Fe-As layer with the coordinate system. 
Note that the center As atom is above the Fe-plane, while the corner As atoms are below the 
Fe-plane. 



structures [9|,|38| using the quasi-atomic minimal basis-set orbitals (QUAMBOs) 39N42| . The 
main idea of the QUAMBO approach is to recover a set of local quasi-atomic orbitals by 
performing an inverse unitary transformation for the "bonding" and "anti-bonding" states. 
In practice, the "bonding states" are some occupied bands or low energy bands which are 
intended to be preserved in the new local orbital representation. The "anti-bonding" states 
are constructed from the orthogonal subspace of the "bonding" subspace and are optimized 
such that resultant QUAMBOs have maximal similarity with the free atomic orbitals. The 
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TABLE I: The hopping parameters between the orbitals of Fel/Mnl and those of its nearest 
neighbour Fe2/Mn2, next nearest neighbour Fe3/Mn3 and the center As atom. The entries of 
A1/B1(A2/B2) mean that Al is the hopping parameter in the Fe d-model and Bl in the Mn 
d-model. A2 and B2 are the corresponding values in the ddp-models. 



orthogonal subspace is spanned by the wave functions in some energy window. It becomes 
complete in the infinite energy or band limit, where succinct closed form for QUAMBOs can 
still be obtainedj41|. In principle, different choices of the "bonding" states and the energy 
window for the "anti-bonding" states will yield different set of QUAMBOs and tight-binding 
parameters, which is reasonable and acceptable since the downfolding of the band structure 
should not be unique. We find that the ddp model will give too many Fe/Mn d-electrons if it 
is obtained by treating all the occupied states and virtual states up to 2eV above Fermi level 
as the "bonding" states in the infinite band limit: 7.3 d-electrons per Fe atom for BaFe2As2 



(close to the analysis of Ref. 43|) and 6.2 for BaMn2As2, while the nominal value would be 
6 and 5, respectively. Such large deviation of the correlated electron number clearly do not 
properly reflect the correct physics of these materials. This problem does not exist for the 
Fe/Mn d-model where the electron filling is fixed by assumption. In order to have a general 
tight-binding model suitable for the correlated local orbital-based approaches, we choose a 
finite energy window around the Fermi level for the generation of the "anti-bonding" states 
and another smaller energy window to control the Bloch bands which will contribute to the 
Fe/Mn d orbitals. We can then construct tight-binding models with reasonable number of 
correlated d-electrons. As a result, we have 6.3 Fe 3d-electrons and 5.3 Mn 3d-electrons for 
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FIG. 4: (Color online) The square of the orbital renornialization factor z^^ as a function of Hund's 
rule coupling J parameter at Hubbard U=2eV, 3eV and 4eV in the Fe d-model (upper panels) for 
BaFe2As2 and Mn d-model for BaMn2As2 (lower pannels) in the paramagnetic state. 



the ddp-model. Those are the initial occupancies that are further reduced due to the self- 
consistent determination of the occupancies within the Gutzwiller approach and the double 
counting corrections, reaching values close to 6 Fe 3d-electrons and 5 Mn 3d-electrons. 

Figure [1] and [2] show the bare band structures ([/ = J = 0) of the d-model and the 
ddp-model for BaFe2As2 and BaMn2As2, respectively. The LDA band structures have also 
been plotted for comparison. One can see that the band structure from the ddp-model 
agrees very well with the LDA result in the low energy window for both systems. The 
Fe/Mn d-model, however, fails to reproduce one low energy unoccupied band near F-point 
which is mainly contributed by the Ba 5d orbitals. The total number of electrons is 6(5) 
per Fe(Mn) atom for the d-model and 12(11) per Fe(Mn) atom for the ddp-model. The 
tight-binding parameters between one Fe/Mn atom and its nearest Fe/Mn atom, second 
the nearest Fe/Mn atom and the nearest As atom are listed in Table [H with the atomic 
configuration and the coordinate system shown in Figj3l The hopping parameters of the 
BaFe2As2 system are overall larger in magnitude than those of BaMn2As2 system, which is 
mainly owing to the fact that BaMn2As2 has a much larger volume than BaFe2As2, although 
the ionic radii of Fe and Mn are very close. The direct Fe-Fe/Mn-Mn hopping parameters 
become smaller in the ddp-model, as expected. The Fe d-model has been compared with 



the one downfolded using maximally localized Wannier function (MLWF) 4^, |45|. Although 



the MLWF approach can achieve better fitting after some fine tuning 461], we find that the 
results reported here are not sensitive to this detail. 
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FIG. 5: (Color online) Same as Fig. [Hbut based on ddp-model. 



IV. COMPARATIVE STUDY OF THE D AND DDP MODEL FOR PM PHASE 

Figure H] and [5] show the variation of the z^-factors for the correlated d-orbitals with 
increasing Hund's rule coupling J and Hubbard U fixed at 2eV, 3eV and 4eV in the param- 
agnetic state of the Fe/Mn d- model and the ddp-model, respectively. On the small J side, 
the 2;^-factors for both models of the compounds stay almost the same. With increasing J, 
the z^-factors start to drop rapidly but remain finite for the two models of BaFe2As2. In 
contrast, the 2;^-factors exhibit a sharp decrease to zero for both models of BaMn2As2, i.e., 
Mott localization for all the orbitals sets in beyond a threshold value of J . This indicates 
that BaMn2As2 is a more correlated system than BaFe2As2. By increasing U from 2eV to 
4eV, the transition of the z^-curves occurs at somewhat smaller J-values for all the model 
calculations. The results show a much stronger dependence on the Hund's rule coupling J 
rather than Hubbard f/, which is typically the case in correlated multiband systems. Com- 
paring the Fe/Mn d- model and the ddp-model, we observe a systematic shift of the 2;^-curves 
to larger J value with the inclusion of the Ba 5d and As 4p orbitals to the Fe/Mn d- model. 
Thus, a larger value of J is needed in a ddp-model to obtain the similar z^-factors as in the 
Fe/Mn d- model. This can be attributed to the screening effect of additional electrons in the 
ddp-model as caused the hybridization of the Fe/Mn 3d orbitals with other, less correlated 
degrees of freedom. 

The general feature of the local d-orbital occupations in the model calculations as shown in 
Fig. [6] and [7] is consistent with that of the z^-factors described above. While the Hubbard U 
parameter tends to polarize the orbitals, the Hund's rule coupling J favors the equalization 
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FIG. 6: (Color online) Orbital occupations for the Mn/Fe d-model with settings same as Fig. HI 
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FIG. 7: (Color online) Local Fe/Mn 3d rbital occupations for the ddp-model with settings same 
as Fig. H 



of the local orbital occupations. All local d-orbitals become half occupied as the system 
approaches the Mott transition for the models of BaMn2Fe2 which has 5 electrons in each 
5 Mn 3d orbitals. For the BaFe2As2 system, there is 6 electrons in each 5 Fe 3d orbitals. 



Hence we expect an orbital selective Mott transition would occur with large U and J|47l. 148 1. 
We did perform model calculations at even larger U and find an orbitally selective transition 
where the z factors of all orbitals except the ^^.j^-orbital vanish, while z^y = 1 as the orbital 
becomes completely filled. This behavior requires however unphysically large values of U. 
Indeed some signature can be identified as one can see that n^o 2 becomes closest to 0.5 and 
z^.2_ 2 drops most rapidly with increasing U and J for the Fe d-model. When comparing the 
results of Mn/Fe d-model and the ddp-model, a similar systematic shift of the local orbital 
occupation curves to the large J side is observed due to the screening and hybridization effect. 
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FIG. 8: (Color online) The square of the orbital renormalization factor z^^ as a function of Hund's 
rule coupling J parameter for the ddp-model of BaMn2As2 (solid symbol) and that of BaFe2As2 
(open symbol) with total number of electrons fixed to 11 per Fe/Mn atom and U=3eV. 



We conclude that BaMn2Fe2 is approaching a Mott insulator with large localized moment. 
Combining our calculations in the magnetically ordered state and in the paramagnetic state 
we obtain J ~ 0.2 — 0.3eV to yields the experimentally observed ordered moment, values that 
are not yet in the Mott insulating regime, yet they are rather close. In contrast, BaFe2As2 
is significantly away from Mott localization, yet there are clearly visible polarizations of 
the orbital populations that demonstrate that the system is in an intermediate regime. To 
demonstrate that this is indeed the case we determined the quasiparticle weight of a system 
with the tight binding parameters (ddp-model) of BaFe2As2, yet with a total electron count 
of 11 per Fe atom as appropriate for BaMn2As2 (see Fig. [8]). Now a Mott transition similar 
to the actual BaMn2As2 calculation occurs, albeit at a somewhat larger value of J. 



V. EFFECTS OF CHARGE SELF-CONSISTENCY 



The variation of the 2;^-factors and occupations of the correlated Fe/Mn 3d orbitals with 
U and J has further been investi gate d by the LDA+Gutzwiller method 30|, |3l|, which is built 
upon the LDA+U approach 49|, l50|]. We use the same correlated orbitals obtained via the 
QUAMBO procedure as in the Fe/Mn d- models for the calculations and achieve the electron 
density self-consistency. As shown in Fig. [9l the z^-curves exhibit a similar trend but vary in 
a way much slower than that in the Fe/Mn d- models, which can again be attributed to the 
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FIG. 9: (Color online) Same as Fig. Hlbut based on LDA+G calculations. 
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FIG. 10: (Color online) Local Fe/Mn 3d orbital occupations based on LDA+G calculations with 
settings same as Fig. [H 



screening effect of the other electrons in the systems. Our results on the BaFe2As2 system are 
in good agreement with those reported in Ref. 22|. Overall, the 2;^-factors for BaMn2As2 
are systematically smaller than those for BaFe2As2. Note that the calculated screened 
U(J) is 3.6(0.76)eV for Fe and 3.2(0.70)eV for Mn based on the constrained random phase 
approximation and the maximally localized Wannier function 5l[, it is fairly reasonable to 
assume that the screened interaction parameters in the models investigated here are also 
very close. Therefore BaMn2As2 is a more correlated material than BaFe2As2. 



VI. COMPARATIVE STUDY OF THE D AND DDP MODEL FOR AFM PHASE 

The ground state AFM phase has also been studied for the models of the two systems. 
Figure [H] shows the behavior of the magnetic moment as a function of J with U fixed at 
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FIG. 11: (Color online) Magnetic moment vs J at U=2eV (squares), 3eV (circles) and 4eV (trian- 
gles) in the Fe d-model (upper panel) and the ddp-model (lower panel) for BaFe2As2. The solid 
symbols indicate the results using Gutzwiller approximation, and open symbols using Hartree-Fock 
approximation. 



2eV, 3eV and 4eV in the stripe- type AFM ground state for BaFe2As2. In contrast to the 
Fe d-model including complete local interactions for LaOFeAs in which the system exhibits 
a sharp transition from PM state to AFM state with large magnetic moment using the 
Gutzwiller approximation 2J] , the models with only density- density type interactions for 
BaFe2As2 describe the PM-AFM transition fairly well as the magnetic moment increases 
rapidly yet smoothly when J rises as indicated by the solid symbols in the upper panel. 
The magnetic moment varies from O.lAfis to l.SSfiB in the range of J = 0.22 — 0.32eV 
with U = 3eV. For reference the experimental magnetic moment is 0.87/1^ for BaFe^Aso j52|. 
With the inclusion of the Ba 5d and As 4p orbitals in the ddp-model, the magnetic moment 
increases slower as seen in the lower panel. The high sensitivity of the magnetic moment 
to Hund's rule coupling J is consistent with the DFT+DMFT calculation results|2o|. For 
comparison. Fig. [H] also shows the magnetic moment of the AFM phase in the Fe d- 
model with Hartree-Fock approximation by the open symbols in the upper panel. Within 
Hartree-Fock approximation, the PM-AFM transition takes place at much smaller J. The 
magnetic moment exhibits quite smooth variation with increasing J at U=2eV. However, an 
abrupt transition from PM to AFM state with large magnetic moment occurs at U=3eV. In 
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FIG. 12: (Color online) Magnetic moment vs J curves for BaMn2As2 with setting as as Fig. [TTl 



comparison, the models of BaMn2As2 yield a very sharp transition from PM to the G-type 
AFM ground state with magnetic moment M ^ 3/2^ as shown in Fig. [121 However, it is 
still reasonable since the experimental magnetic moment is 3.88/xs 9|. The z^-factors for the 
models of BaFe2As2 in the AFM ground state are always or order unity (> 0.8). 



VII. CONCLUSION 



In summary, we have studied the Fe/Mn d-model and the ddp-model on the electronic 
and magnetic properties of the BaFe2As2 and BaMn2As2 systems. The renormalization 
factors of the correlated Mn 3d orbitals are found to be systematically smaller than those of 
Fe 3d orbitals, implying that the electron correlation for BaMn2As2 are much more efficient 
to cause Mott localization physics compared to BaFe2As2. Ultimately this is due to the fact 
that Fe, with its close to six electrons, must undergo an orbital selective Mott transition while 
the odd number of electrons in case of Mn allow for a more ordinary Mott transition. While 
the strength of the interactions are not sufficient for Mott localization in either system the 
Mn-based material is significantly closer to localization. The LDA+Gutzwiller results on the 
paramagnetic phase with electron density self-consistency also confirm the conclusion. The 
variation of the magnetic moment in the AFM ground state of the two compounds seems in 
accordance with the experimental results: a smooth increase of the magnetic moment with 
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rising J for the stripe-AFM ground state of BaFe2As2 with small experimental magnetic 
moment and an abrupt jump of the magnetic moment from to about S/ib for the G-AFM 
ground state of BaFe2As2 with large experimental magnetic moment. We also checked that 
the two different ordered states are indeed energetically lower in the respective systems. 
The Gutzwiller approximation, however, is not able to provide a quantitatively consistent 
description of the AFM ground state for the iron pnictide systems with both correct magnetic 
moment and band renormalization factor for the same parameters. Nevertheless it is a 
comparatively easy to implement and powerful tool that allows for a first analysis of the role 
of magnetic correlations, the nature of magnetically ordered states in strongly correlated 
multi-orbital materials. 
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